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ABSTRACT 

We have detected weak soft X-ray emission from the Pnlsar Wind Nebnla trailing the high 
velocity star PSR 2224-1-65 (the ‘Gnitar Nebnla’). This X-ray flnx gives evidence of 7 ~ 10^ 
eV particles in the pnlsar wind and constrains the properties of the post-shock flow. The 
X-ray emission is most easily nnderstood if the shocked pnlsar wind is partly conhned in the 


nebula and if magnetic helds in this zone can 

1. Introduction 

Identifying the spindown products of rotation-powered 
pulsars is a long standing problem in astrophysics. Recent 
observations of high energy emission from young pulsars 
(eg. Fierro 1996) and improvement in emission models 
(eg. Romani 1996) show that up to a few xl0% of the 
spindown power can be ascribed to observed photons. In 
most cases, the remaining power in particles and unob¬ 
served radiation still needs to be identified. The wind of 
energetic particles and fields emitted from spin-powered 
pulsars is most easily observed when it is confined. Such 
post-shock, or ‘plerionic’ emission has long been associated 
with sources such as the Crab and Vela pulsars and with 
similar filled-center remnants such as 3C58, presumably 
containing young pulsars. 

The discovery of bow shock nebulae around a number of 
pulsars has provided a further opportunity to study pulsar 
winds. These nebulae generally occur when the pulsar is 
either long-lived, with a low magnetic field B and a short 
spin period, or when it has a large space velocity; these 
properties allow escape from the pulsar’s parent supernova 
remnant while it retains substantial spindown power. The 
pulsar wind then shocks against the ambient ISM. These 
objects are most prominent as Ha nebulae, first discov¬ 
ered around the millisecond pulsar PSR 1957-1-20 (Kulka- 
rni and Hester 1988, see also Aldcroft, Romani and Cordes 
1992). A similar bow shock is seen around the 5 ms pul¬ 
sar PSR J0437-4715 (Bell, Bailes and Bessel 1993) and 
Cordes, Romani and Lundgren (1993, CRL) have detected 
an Ha nebula around the high velocity, slow spin pulsar 
PSR 2224-f65. 

High energy observations of these pulsar wind nebulae 
(PWNs) can place key constraints on the pulsar particle 
output. Models of the Crab nebula (Kennel & Coroniti 
1984), for example indicate that the bulk of the spindown 
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grow to near equipartition values. 


power is carried off in an wind with 7 ~ 10®'^. Arons 
and Tavani (1993) described the expected high energy 
emission from wind shocks in the PSR 1957-1-20 system, 
focusing on a wind model including a massive baryonic 
component. Using ROSAT Position Sensitive Proportional 
Counter (PSPC) observations of PSR 1957-1-20, Kulkarni 
et al. (1992) were able to extract general constraints on 
the content of the pulsar wind. A trail of soft X-ray 
emission behind the nearby pulsar PSR1929-I-10, discov¬ 
ered in ROSAT PSPC data by D. Helfand and described 
by Wang, Li and Begelman (1993), is also a likely product 
of pulsar wind-ISM interaction. Finally the interaction of 
pulsar winds with companion stellar winds can also pro¬ 
duce a bow shock geometry. For example, PSR 1259-63 in 
an eccentric orbit around a Be companion provides the op¬ 
portunity of studying the pulsar wind interactions under 
changing conditions (Arons and Tavani 1994). ASCA ob¬ 
servations of non-thermal X-rays (see Hirayama 1996, and 
references therein) have in particular provided important 
constraints on the particle energy and luminosity. 

We describe here ROSAT High Resolution Imager (HRI) 
observations of the PWN trailing PSR 2224-1-65, the ‘Gui¬ 
tar nebula’. This is followed by a description of the bow 
shock flow geometry and a discussion of the origin of the 
observed X-ray flux. The Guitar nebula is particularly in¬ 
teresting because the pulsar itself is relatively normal, but 
the high velocity and peculiar shock geometry of the sys¬ 
tem provide opportunities to probe the pulsar wind inter¬ 
actions under conditions quite different than those found 
in other PWN. In particular, the small bow shock scale 
allows soft X-ray constraints on electron energies similar 
to those inferred for the Grab pulsar. 
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2. ROSAT HRI flux 

We proposed ROSAT HRI observations of PSR 2224+65 
to detect emission in the shocked wind and ISM and to 
constrain the physics of the post-shock flow. Exposures 
were obtained in 7 observation intervals during July 13 - 
25, 1994. With 36,894s of dead-time corrected exposure 
centered on PSR 2224+65 we were able to detect coro¬ 
nal emission from ~ 10 bright stars in the field. Cross 
correlation with optical frames registered to deep Ha im¬ 
ages showing the Guitar nebula allowed us to ‘boresight’ 
the HRI data. Three Hubble Space Telescope (HST) guide 
stars (CSC 472-144, -188, -232) were detected with stan¬ 
dard processing in the center of the HRI frame and were 
used for final alignment of the X-ray and optical images 
to < 1". This allowed a very sensitive search for emis¬ 
sion from PSR 2224+65 and the Guitar Nebula. No X-ray 
emission was detected from PSR 2224+65 itself, with the 
2 counts detected within 4" of the pulsar position entirely 
consistent with the 2.02 counts expected from the mean 
background rate. However an X-ray excess was appar¬ 
ent in the vicinity of the Ha nebula. We have defined 
an aperture of half width 4" following the outline of the 
edge-brightened Guitar shape. If the X-ray emission is 
confined to the limb, then this contour will contain 90% of 
the source counts; if the bubble interior has X-ray emis¬ 
sion, the omitted flux will be slightly larger. This contour 
contains 31 counts at the Guitar position in excess of the 
58 counts expected from the mean background rate (Figure 
1). This is formally a 4.1 ct detection. To check the signifi¬ 
cance, we summed counts in our nebula-shaped aperture at 
various independent positions in the HRI image, after ap¬ 
plying vignetting corrections and excluding regions where 
significant point sources (generally identified stars) were 
found in the standard PROS processing. The histogram 
of the resulting excess counts shows that the true Guitar 
position has the largest excess. Estimating the envelope 
of the count fluctuations again gives a significance slightly 
in excess of 4 ct. The Guitar nebula thus emits soft X-rays 
with a count rate of 9.6 x 10“‘*counts/s in the ROSAT 
band. The flux seems to follow the limb of the Guitar, 
concentrating in the forward facing parts of the nebula, 
but starting significantly behind the pulsar itself. 

Some crude spectral information is available from the 
PHA energy channels in the HRI detector. We find that 
the excess Guitar counts are soft, with a median PHA 
channel of 6. The distribution is only slightly harder than 
that of the mean background counts, but also contains an 
excess of counts in PHA channels 0-1, suggesting UV flux 
near the source. This ‘spectrum’ is clearly different from 
some of the stronger sources identified with bright stars in 
the optical; these have median PHA channels of 8-9 and 
very few counts below channel 5. Unfortunately, with the 
very low count rate and lack of calibration, quantitative 
results cannot be extracted from this hardness informa¬ 
tion. 

The implied X-ray luminosity is sensitive to the par- 



Fig. 1.— Above: Ha image of the Guitar Nebula, showing 
our shaped aperture. Below: Aperture over central portion 
of the smoothed HRI count map. 


ent spectrum and the absorption. From the pulsar disper¬ 
sion measure 35.3cm“^/pc and an assumed ne/nn = 0.1 
we infer a column density Nh ~ 10^^cm“^. This as¬ 
sumed UelnH is very similar to that inferred for PSR 
2334+61 (Becker, Brazier & Triimper 1996) at a similar 
direction and DM. The Balmer decrement in the Gui¬ 
tar nebula itself suggests Ay = 1.65, corresponding to 
Nh = 3 X 10^^cm“^. However, because the Balmer ratios 
in non-radiative shocks are uncertain (Chevalier, Kirsh- 
ner and Raymond 1980; Raymond 1991), we adopt here 
10^^ for the absorbing column. With this Nh and a pho¬ 
ton number index of a = 2 we infer an X-ray flux of 
1.8x 10“^^erg/cm^/s (O.Ol-lOkeV); for A// = 3xl0^^cm“^ 
the flux is 50% higher. This represents 1.7 x 10“^(i^p^//45 
of the spindown luminosity. 

3. Geometry of the Guitar Nebula 

As described in CRL, the guitar nebula is an Ha- 
emitting bow shock trailing the slow spin (P=0.68s) inter¬ 
mediate age (tc = 1.1 X 10®y) pulsar PSR 2224+65. With 
an interferometrically determined proper motion giving a 
space velocity of Vp = 8.4 x lO^dkpc/cosi cm/s (with i the 
angle between the velocity and the plane of the sky) and a 
dispersion measure distance estimate of dkpc = 1-95 kilo- 
parsecs (Taylor & Cordes 1993) this object has arguably 
the highest space velocity of any known star. 
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Fig. 2.— Histogram of excess counts in non-overlapping 
Guitar-shaped apertures in the HRI image. The observed 
counts at the Guitar position are in the shaded bin. 

With this velocity the standoff of the pulsar wind bow 
shock is especially small for this nebula. The pulsar spin- 
down luminosity E = 1.3 x 10^^745 erg/s is believed to 
power a relativistic wind which shocks against the ISM 
of density p = 2.2 x g/cm^, giving a momentum 

balance standoff of the contact discontinuity at the bow 
shock nose of 

^0 = ( . ^ ^ )^/^cosz/(ikpcCm. (1) 

The Balmer emission of the pulsar wind nebulae is caused 
by neutral ISM atoms drifting into the shocked gas and 
suffering collisional excitation or charge exchange into ex¬ 
cited states of the post-shock ions. Accordingly Ha arises 
in a thin sheet behind the ISM shock, upstream from the 
contact discontinuity (Ghevalier, Kirshner and Raymond 
1980). A detailed model for the bow shock flow is not avail¬ 
able, but Aldcroft, Romani & Cordes (1992) estimated 
that the Ha apex lies at ~ 1.3ro. This size was unre¬ 
solved in ground based images of the Guitar. HST Wide 
Field Planetary Camera (WFPG2) Ha images of the neb¬ 
ula, however, show that the apex has the classic bow shock 
shape displayed by the PSR 1957-1-32 and PSR J0437-4715 
nebulae. These data are described elsewhere (Cordes et al. 
in preparation). Here we note only that while the pulsar, 
and hence the standoff angle, are not directly observed in 
the HST data, a fit to the bright Ha ridge of the bow shock 
allows us to measure the size of the bow shock and infer 
tq. Wilkin (1996) has derived convenient analytic expres¬ 
sions for the geometry of a thin momentum-balance bow 
shock. Assuming that the shock standoff at the apex scales 
proportionally for the rest of the bow wave, the scaled ver¬ 
sion of the momentum balance contact discontinuity curve 
gives the bow shock shape: 

OhcW = l-3r((())/d = 1.3ro csc(()\/3(l — ^cot (l))/d (2) 

where the pulsar is at the origin and position angle </> is 
measured from the pulsar velocity. The fit determines the 


origin and the overall scale, giving an inferred Ha standoff 
angle as OnaiO) = 0.038" ± 0.005". This corresponds to 
a contact discontinuity scale rg = 4.3 x lO^^dkpc/cosicm 
which, with ( 1 ) gives an estimate of the pulsar distance as 

dkpc = l.lcost (/45/nR)^/'‘- 

Ground-based Palomar 200" images (Figure 1) show 
that the shape of the Ha nebula downstream from the 
HST resolved apex can be described as a ~ 3" x 20" tube 
(the ‘neck’ of the Guitar) connecting to a complex bub¬ 
ble that flares gradually to a diameter of ~ 30" some 70" 
behind the pulsar. The waist in this bubble (the Guitar 
‘body’) may be due to ISM density variations. The bub¬ 
ble also does not appear fully closed at the end opposite 
the pulsar. With the observed proper motion, the pulsar 
traverses the Guitar in ~ 430y. CRL give the Ha flux of 
the head of the nebula as 1.1 x 10“^photons/cm^/s. This 
line flux implies that neutral H atoms enter the nebula 
head at 1.3 x 10^^d|p^atoms/s. With the observed radius 
of ~ 1.5", an ISM neutral fraction X and the measured 
pulsar proper motion, the nebula head should sweep up 
neutral H at 1.3 x lO^^ni/dlp^/Aatoms/s. Thus the ob¬ 
served line flux is in good accord with the HST constraints 
on the bow shock geometry for nndkpc/X ~ 1 . 

The variation in nebula cross section can be understood 
from pressure balance arguments. Behind the pulsar the 
outflowing wind drives the bow shock tail into the ISM; 
this stalls where the wind ram pressure balances the ISM 
static pressure ~ 10 “^^P_i 2 g/cm/s^. Well behind the pul¬ 
sar the wind impacts at a shallow angle; a characteristic 
position angle of ~ 150° gives a transverse pressure ~ 
Esm^{TT — 4>)/{4TT0‘j_d^c) ~ 2xlO“^^/45(6>X'^fcpc)~^g/cm/s^ 
in a tube of radius d" arcseconds. Thus we expect the 
ram pressure driven ‘neck’ of the Guitar to stall at a 
width 6± ~ 1 . 4 "(/ 45 /P_i 2 )^/^/(ifcpc- Further downstream 
it is clear that the shocked pulsar wind is partly con¬ 
fined in the body of the Guitar. This confinement may 
be aided by the increased ram pressure associated with 
expansion into a dense region toward the rear of the Gui¬ 
tar, delimited by the horizontal filament in Figure 1. Here 
the isotropic pressure of the shocked relativistic plasma 
Prei = u/3 comes into play. If the wind emitted dur¬ 
ing the time the pulsar crosses the Guitar is fully con¬ 
fined in a tube of diameter rj_ and there is no cooling, 
the pressure is E/(Sirr'^Vp). In practice the back of the 
nebulae does not appear closed, so perhaps ~ 0.1 of this 
pressure obtains; even with this pressure the bubble stalls 
with a characteristic radius 9± « 26 "(/ 45 /P_i 2 )^/^dfcp/^. 
We thus have a consistent picture of the Guitar bubble 
as a classic bow shock whose ram pressure stalls against 
the static ISM pressure behind the pulsar; the backflow 
in this neck should be mildly relativistic. Further down¬ 
stream, the shocked pulsar wind stalls and becomes partly 
confined. The shocked wind expands against the ISM in¬ 
flating the bubble to the observed ~ 30" diameter. 
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4. Implications and Bow Shock Flow 

We turn now to the interpretation of the lack of X- 
ray flux from PSR 2224+65 and to the X-ray flux from 
the Guitar nebula body. The non-detection of the pul¬ 
sar is not surprising. At Tc = 10®y the expected surface 
temperature of ~ 3 x lO^K from standard cooling models 
predicts a luminosity of ~ 6 x 10 ^° erg/s and no signifi¬ 
cant HRI flux. Pulsars of similar age have been detected 
by ROSAT, eg. PSR 1929+10 (Yancopoulos et al. 1994) 
and PSR 0823+26 (Sun et al 1993), but the emission from 
these nearby objects seems to be from a small heated po¬ 
lar cap; scaled to PSR 2224+65’s distance the polar cap 
flux of PSR 1929+10 would contribute < 1 count in our 
exposure. 

The two plausible sources of the nebular X-ray emission 
are the shocked pulsar wind and the postshock ISM. In 
contrast to the bow shocks around the slower pulsars PSR 
1957+20 and PSR J0437-4715, the large pulsar velocity 
here gives high T in the shocked ISM gas; at the apex of the 
bow wave, the postshock ISM temperature is 3 x lO^d^p^K. 
However, the small working surface resulting from the 
small ro means that the shock luminosity near the Guitar 
nose is only ~ T:rlpv‘pl2 ~ 8 x lO^^n^^d^p^erg/s. More¬ 
over, for a typical ISM abundances and cooling rate A, the 
characteristic cooling time 3fcr/(8n_f/A) ~ 8 x lO^^nj^^s 
is quite long, so little thermal emission is expected be¬ 
fore adiabatic losses in the post-shock flow become signif¬ 
icant. The absence of efficient cooling is supported by the 
Balmer-dominated nebular spectrum at the apex (GRL) 
which indicates that the shock is non-radiative (Gheva- 
lier, Kirshner and Raymond 1980). Along the body of the 
Guitar, the continued expansion of the partly confined pul¬ 
sar wind drives shock waves into the surrounding medium 
at w ~ 1.5 X lO^dfcpcCm/s. The implied shock luminosity 
along the Guitar body is Lst ~ 2 x 10^^d^p^n//erg/s, with 
a postshock temperature of ~ 4 x lO^d^p^K. Although the 
cooling time for d^pc ~ 1 is comparable to the flow time 
through the nebula (so that we might expect some thermal 
emission from the tail of the Guitar), the effective temper¬ 
ature is too low to produce the observed HRI count rate. 
Even if the distance to PSR 2224+65 is as large as 2kpc, 
prospects for thermal emission from the nebula are not 
much better. With T^h = 1-6 x 10®K our observed HRI 
counts imply a thermal luminosity at 2 kpc of 8 x lO^^erg/s 
for Nh = 10^^cm“^ and a Raymond-Smith thermal spec¬ 
trum. However at these temperatures the cooling time is 
~ 3 X so that radiation in the postshock flow 

is inefficient. Thus unless dkpc and nn are large, Nh is 
small and the cooling rate is significantly enhanced, ther¬ 
mal emission cannot account for the observed HRI flux. 

A more promising source is post-shock synchrotron 
cooling of the relativistic pulsar wind. In the ram-pressure 
balance region of the bow wave, the energy dissipated in 
this shocked wind is larger than that in the shocked ISM 
by c/vp ~ 300. Observations of the Grab nebula plerion 
suggest that the pulsar wind is dominated by energetic 


with 7 = 10^77 , 77 ~ 0.3. These electrons radiate in the 
postshock field of the pulsar wind, Bg = 2x 10“^(ro/r)G. 
Kulkarni et al. 1992 were able to use the X-ray flux from 
the millisecond pulsar PSR 1957+20 to limit the fraction 
of the pulsar spindown in with 77 ~ 0.01 and 77 ~ 10 , 
assuming synchrotron emission at the companion shock 
and at the bow shock, respectively. With it’s larger field 
and high velocity, PSR 2224+65’s soft X-ray flux is most 
sensitive to intermediate electron energies, similar to those 
in the Grab Nebula. Including the r~^ dependence of the 
pulsar wind B we find that the synchrotron spectrum in 
the shocked wind has a peak energy Ec ~ 3 . 5 (ro/r) 77 keV, 
so that closer than rmax ~ the charges will radiate 

in the ROSAT band. 

We expect the shocked pulsar wind to have a mildly 
relativistic backflow in the bow shock, with a charac¬ 
teristic flow time at a distance r from the pulsar of 
Tfi ~ 3r/c. This is shorter than the synchrotron cool¬ 
ing time under these conditions, so that the efficiency is 
small rj ~ Tfilrgy ~ 3 x 10 “^(ro/r) 77 . When the peak 
energy is well above the ROSAT threshold the observed 
flux will decrease somewhat; since the synchrotron emis- 
sivity at energy E scales as ~ A^/^exp(—A/Ac), the frac¬ 
tion of the radiation in the ROSAT band has the weak 
dependence (0.3keV/i?c)^/^ ~ 0 . 4 (r/ro)^/^ 7 '^^^^. Thus 
if a fraction of the spindown energy is carried off by 
charges of Lorentz factor 7 , then including the cooling ef¬ 
ficiency and the fraction of the spectrum in the ROSAT 
band, we get an observed X-ray luminosity from the wind 
shock within position angle (p of 


L, 


10-V7-E(1 - cos 4 >) 


IdL 

[r{(j))\ 


2/3 


1/3 

77 


where r{4)) is determined from (2). Thus the emission 
should be concentrated toward the bow shock apex. Be¬ 
cause the wind cannot cool effectively before the relativis¬ 
tic backflow carries it towards the body of the Guitar, the 
total luminosity is small. With an a = 2 spectrum and 
Nh = 10 ^^cm“^ at Ikpc, we expect 2 counts from the 
shocked pulsar wind at the apex of the flow. 

Above we have assumed that the field in the pulsar wind 
itself dominates B in the post-shock region and that the 
flow is not confined. However, it appears that the pulsar 
wind is partly confined in the nebula. In addition tur¬ 
bulence in the shocked pulsar wind may increase B up to 
equipartition values. At the shock apex Bgq Bg, so that 
our previous estimates hold. However, ~ 30ro (~ 1") be¬ 
hind the pulsar the ram pressure associated with the side¬ 
ways expansion of the guitar body at 1.5 x IttddkpcCm/s be¬ 
gins to exceed B"^/Att in the shocked pulsar wind. Thus as 
the neck of the Guitar merges into the body we expect that 
an equipartition field can reach Bgq ~ 8 x IQ-^dkpcnJp'^G. 
The characteristic energy in this field Ec ~ 0.1577keV is 
large enough to contribute significant ROSAT emission. 
More importantly, if the shocked wind is partly confined in 
this region, we expect t/; ~ ~ 5 x 10®s. For dkpc = 1 

we get a synchrotron cooling time Tgy ~ 9 x 10 ® 7 ^^s and 
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about 30% of the shock energy radiated in the ROSAT 
band for 77 = 1. This gives an observed X-ray flux 
0.3EfjTfi/Tsy Thus to produce the observed X-rays, we 
require ~ 0.1 in 77 > 1 particles. With the increased 
Reg inferred at 2 kpc, cooling is even more efficient and we 
again infer « 0 . 1 . 

5. Conclusions 

The X-ray emission observed from the bow shock cre¬ 
ated by PSR 2224-1-65 represents ~ 0 . 02 (i^p^/J 45 of it’s 
spindown power. Although the pulsar’s large velocity en¬ 
sures that gas is heated to X-ray emitting temperatures 
near the bow wave, the cooling in the post-shock ISM flow 
is inefficient. Little of the observed flux can be contributed 
by the shocked ISM. The location of the X-ray emission 
significantly downstream from the pulsar suggests that the 
most likely source of the flux is the shocked pulsar wind. 
However, to account for the observed flux we require this 
wind to be partly confined in the Guitar body and that 
turbulence in the shocked wind flow can amplify the mag¬ 
netic field to near equipartition vales. If these conditions 
obtain, the observed flux implies that > 0.1 of the pul¬ 
sar spindown power is being emitted in with 7 ~ 10 ^. 
Our data then support the idea that a pair plasma with 
10 ® - 10 ^ is a dominant component of the pulsar 
wind for normal non-recycled pulsars. Such a wind is in 
reasonable accord with models of pulsar magnetosphere 
acceleration, but the pairs may also be accelerated in the 
post-shock region itself. In contrast the conditions inferred 
for the soft X-ray emission from PSR 1929-|-10’s shock by 
Wang, et al. suggest that 7 > 10® pairs are responsible for 
the PSPC photons; it is perhaps significant that the flux 
from these pairs represents only 10“"^ of PSR 1929-|-10’s 
spindown luminosity. 

A more complete description of thermal and non-thermal 
emission in the Guitar nebula will require detailed hydro- 
dynamic models of the backflow in the postshock region. 
For example, Raga, Cabrit & Canute (1995) find that there 
is mixing of the shocked fluids in a bow shock flow. While 
the presence of the relativistic pulsar wind means that 
their sums are not directly applicable to PSR 2224-1-65, 
such mixing could affect the estimates of cooling time in 
the postshock flow. Such modeling coupled with future 
high throughput X-ray observations could give important 
constraints on the pulsar wind in this unique system. 

This work was supported in part by NASA grants NAG 
5-2706 and NAGW-4526 (RWR), and NSF grant AST 92- 
18075 (JMC). 
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